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Abstract 

Helium and hydrogen recombination lines observed in low-metallicity, extragalactic 
Hii regions provide the data used to infer the primordial helium mass fraction, Yp. The 
ionization corrections for unseen neutral helium (or hydrogen) are usually assumed to be 
absent; i.e., the ionization correction factor is taken to be unity (icf = 1). In this paper we 
revisit the question of the icf for Hii regions ionized by clusters of young, hot, metal-poor 
stars. Our key result is that for the Hii regions used in the determination of Yp, there is a 
"reverse" ionization correction: icf < 1. We explore the effect on the icf of more realistic 
inhomogeneous Hii region models and find that for those regions ionized by young stars, 
with "hard" radiation spectra, the icf is reduced further below unity. In Monte Carlos 
using Hii region data from the literature (Izotov and Thuan 1998) we estimate a reduction 
in the pubhshed value of Yp of order 0.003, which is roughly twice as large as the quoted 
statistical error in the Yp determination. 
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1. Introduction 

The primordial abundance of ^He is key to testing the consistency of the standard 
hot big bang model of cosmology and to using primordial nucleosynthesis as a probe 
of cosmology and particle physics (for a recent review, see Olive, Steigman & Walker 
1999). Since stars burn hydrogen to helium in the course of their evolution contaminating 
any primordial helium in the interstellar gas with their debris, it is the availability of 
large numbers of carefully observed, very low-metallicity extragalactic Hii regions which 
has permitted estimates of the primordial helium mass fraction, Yp, whose statistical 
uncertainties are very small, 1% {e.g., see, Olive & Steigman 1995 (OS), Olive, Skillman 
& Steigman 1997 (OSS), Izotov, Thuan & Lipovetsky 1994, 1997 (ITL), Izotov & Thuan 
1998 (IT)). However, along the path from the observational data to the derived abundances, 
contamination by unknown systematic uncertainties may bias the inferred value of Yp. 
Observers have identified many potential sources of systematic error (Davidson & Kinman 
1985; Pagel et al. 1992 (PSTE); Skillman et al. 1994; Peimbert 1996; ITL; IT; Skillman, 
Terlevich & Terlevich 1998) and, where possible, have designed their observing programs to 
minimize and/or to account for them. In a previous paper (Steigman, Viegas & Gruenwald 
1997) we have explored the magnitude of the contribution to systematic error in Yp from 
possible temperature fluctuations in the extragalactic Hii regions, concluding that they 
may have a significant effect on the determination of the primordial abundance of helium, 
comparable to or even greater than the statistical uncertainties. Here we turn our attention 
to another potential source of systematic error: the ionization correction for unseen neutral 
hydrogen and/or helium in Hii regions. 

Following Peimbert & Costero (1969), an empirical method is usually employed to 
derive the abundances from the observed emission-line intensities. The electron density and 
temperature of the gas are obtained from various emission-line intensity ratios and they 
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are then used to calculate the appropriate line emissivities which, along with the observed 
emission line intensities, provide the fractional abundances of the various ions which are 
finally combined to obtain the element abundances (see, e.g., Osterbrook 1989). Since not 
all the ions present in the gas actually produce observable lines, an ionization correction 
must be applied to account for the "missing" ionization states. Historically, the ionization 
correction factor (icf ) has been derived for H, He, and many of the heavier elements 
from considerations of the ionization potentials or from numerical photoionization models 
(Peimbert and Torres-Peimbert 1977, Stasinska 1980,1982, Mathis 1982, Pena 1986). In 
the quest for Yp a key concern has been that due to its higher ionization potential, unseen 
neutral helium may lurk in those parts of the H ii regions where the hydrogen is still fully 
ionized {icf > 1). As the metaUicity of the stellar population responsible for creating the 
Hii region decreases, the stars providing the ionizing radiation are expected to be hotter, 
with "harder" spectra. If so, the helium and hydrogen ionization zones should be more 
nearly coincident and the ionization correction minimized {i.e., the icf should be closer to 
unity) . 

Some years ago, Pagel et al. (PSTE) proposed a method for estimating the helium 
icf based on the "radiation softness parameter" 77 = (0+/S+)(S''"+/O++) of Vilchez & 
Pagel (1988). Comparing with photoionization models, they concluded that the icf will 
differ negligibly from unity for log r] < 0.9, corresponding to models with effective stellar 
temperatures higher than 37 000 K. When using a sample of Hii regions to derive Yp, 
PSTE impose this rj condition to discard from their data set those H 11 regions for which 
the icf may differ from unity (log r] > 0.9) and they adopt icf = 1 for all those Hii regions 
with log r) < 0.9. In fact, as we shall see, for the sufficiently "hard" radiation field provided 
by a population of young, hot stars, the ionization correction may actually be reversed with 
neutral hydrogen present in those parts of the H 11 regions where the helium is still ionized 
{icf <1) (Stasinska 1980,1982, Pena 1986, Dinerstein & Shields 1986). In either case. 
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a potentially serious hidden assumption is that of homogeneity for the H ii regions (see, 
e.g., Dinerstein & Shields 1986). Recent HST imaging reveals that these regions are far 
from homogeneous, showing many different features such as pillars, globules and even voids. 
As a result the true icf may differ significantly (on the scale of the statistical errors) from 
the values calculated in simplified models of homogeneous regions, introducing a systematic 
error in the calculation of the primordial abundance of ^He. For example, for Yp 0.24, if 
the icf should differ from unity by 1 - 2%, the change in Yp will be of order 0.002 - 0.004, 
comparable to the statistical uncertainties in Yp suggested by the analyses of OS, OSS, and 
IT. Furthermore, if indeed the icf is less than unity for sufficiently young and metal-poor 
Hii regions, then even the seemingly conservative assumption of PSTE to adopt icf = 1 
may lead to a (metallicity-dependent) bias which systematically overestimates Yp . 

Another common simplification of the photoionization models used to estimate the 
icf is the assumption that the ionizing radiation is provided by stars of a single stellar 
temperature. In contrast, real Hii regions are ionized by clusters of stars of differing masses 
and temperatures whose ionizing radiation spectra differ from that of a single blackbody. 
In the following we adopt the (time-dependent) spectrum of a starburst (Cid-Fernandes et 
al. 1992) appropriate to low-metallicity stars and employ a numerical photoionization code 
(AANGABA; Gruenwald & Viegas 1992) which allows us to account for inhomogeneity. 
The models are described in §2. The icf results are discussed in §3 and are applied to an 
Hii region sample taken from the literature in §4. Our concluding remarks appear in §5. 

2. Photoionization Models 

Many independent 1-D photoionization codes have been developed over the last 30 
years or so. Comparison among several of them for some standard test cases of Hii regions, 
planetary nebulae, and active galactic nuclei reveals very good agreement (Pequignot 
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1986, Ferland et al. 1995). These codes are in wide use analyzing observed emission-line 
spectra. For such codes the usual input parameters are the ionizing radiation spectrum, 
the gas density, and the relative chemical abundances. In most of the codes spherical or 
plane-parallel symmetry is assumed and the diffuse radiation is treated in the outward-only 
approximation. Usually a constant density is assumed, or a specific functional form for 
the variation of density with distance is adopted. However, these simplifying choices may 
not reflect the true structure of real H ii regions as revealed by HST imaging. One way to 
approach the problem of more realistic modeling is to mimic the presence of condensations 
or voids by combining models for several different choices of input parameters. 

In our analysis the AANGABA photoionization code is used to model the gas whose 
distribution is taken to be spherically symmetric. Later, in constructing more general 
inhomogeneous models, this restriction to spherical symmetry is relaxed. Because we are 
specifically interested in the helium/hydrogen icf with the goal of analyzing its effect on 
the helium abundances derived from the low-metallicity H ii regions, a metal-poor chemical 
composition (0.1 solar) is chosen. The stellar cluster ionizing radiation spectrum from 
Cid-Fernandes et al. (1992) is adopted and models are built for two evolutionary phases of 
the star cluster: the initial phase (t = 0), when the spectrum is dominated by the hottest, 
most massive stars, and a later phase (t = 2.5 Myr), when the hot stars have evolved and 
there are fewer He^ ionizing photons. In order to explore a variety of different physical 
conditions, the radiation intensity is characterized by the number of ionizing photons above 
the Lyman hmit, Qi/. Qij is directly related to the mass and IMF of the stellar cluster 
and, along with the gas density, defines the H ii regions. In terms of the popular ionization 
parameter, U, for a given density each of our models corresponds to a fixed value of URj^ 
which is proportional to Qij (Rj is the inner radius of the Hii region). As a result, our 
constant density models should reproduce the icf - rj behaviour found in PSTE (their 
Figure 6). Since the ionization parameter is proportional to Qh/i^, a grid of models at a 
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given density and characterized by the value of the ionization parameter, corresponds to 
different values of Q^^. It is true that those combinations of Qi^, n, and Rj which keep 
U constant lead to the same ionization structure of the gas. However, when looking for 
models to represent the observations of a specific H ii region, ionized by a given number of 
stars, any conclusions should be based on a grid with Qh constant. Thus, in the following 
we show the results of a series of models with different values of U, satisfying the condition 
Qh — constant. First, for each choice of Qh, a homogeneous (H) constant density model is 
constructed. 

The above homogeneous models form the basis for the inhomogeneous models we 
create next: models with spatially limited density enhancements which will mimic density 
condensations (C) and those with density deficits for "valleys" (V). These alternate models 
all start with the same input parameters as the corresponding homogeneous model. We 
set the location of the condensations or valleys by specifying the corresponding optical 
depths at the Lyman limit, r, in the homogeneous model. In our analysis we placed 
the condensations/valleys at the following locations: r = 0.02, 0.03, 0.04, 0.10, 0.40. At 
each of these radii the density (in a spherical shell) is either increased by a factor of 5 
(condensation) or decreased by a factor of 2 (valley), extending over a distance which is 
10% of the thickness of the H ii region for the corresponding homogeneous model. To avoid 
numerical problems, the increase/decrease of the density followed a smooth analytical form. 
Notice however that the size of this transition zone is always much smaller than the size of 
the condensation or the valley. Aside from the condensation/valley, the density throughout 
the rest of the region is the same as in the corresponding homogeneous model. We explored 
many other choices of these parameters but the ones adopted here provide a fair sample of 
the variations of the ionization correction factor (see the next section) in all the models we 
constructed. By constructing inhomogeneous H ii regions using as building blocks a variety 
of the H, C, and V models, we can mimic the physical conditions in realistic Hii regions 
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using the Monte Carlo method. 

Examples of a different sort of composite model are those of Dinerstein & Shields 
(1986), and of Pena (1986) in which they account for the contributions from superposed 
H II regions with high and low stellar temperatures ( "hard" and "soft" radiation fields) . In 
our analysis we can do the same, combining models for H ii regions using the stellar cluster 
ionizing radiation spectra at t = and at t = 2.5 Myr. In the following we pursue both 
kinds of inhomogeneous models. 



3. The Helium Ionization Correction Factor 

In H II regions often only lines from the recombination of He"*" are observed, although 
occasionally those from the recombination of He"'"'^ are also seen. In order to account for 
the possible presence of any unseen He+'^ and He°, and also of H°, we define the icf so that 
the He/H abundance ratio is given by He/H = icf x [n{He'^)/n{H'^)] so that, 

^^^^ W£f!)+^,/[,^^,. 

In addition, following Vilchez & Pagel (1988), the radiation "softness" parameter is defined 
as, 

rj = (n(0+)/n(5+))(n(5++)/n(O++)). (2) 
Both the icf and r) are obtained from the models described in the previous section. 



3.1. Radiation Softness Parameter 

The use of the radiation softness parameter has been challenged by Skillman (1989). 
Using photoionization models characterized by the value of the ionization parameter, he 
found that "/or T^ff > 45,000 K, rj shows a strong dependence on U\ This conclusion could 
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invalidate the use of r] to obtain the hehum ionization correction factor. Skillman's models 
for different H ii regions are characterized by different values of the ionization parameter. 
We note that Skillman uses the Dinnerstein & Shields (1986) ionization parameter which 
is defined as the ratio beween the ionizing photon density and the gas density evaluated 
at the Stromgren radius when the inner radius is fixed at 10% of the Stromgren radius 
{Ri/Rg = 0.10). Each choice of Skillman's ionization parameter corresponds to a different 
Hii region, irradiated by a different immheic of ionizing photons Qh, leading to a different 
ionization structure. It is, therefore, not surprising that when Skillman varies the ionization 
parameter he finds variations in the radiation softness parameter. We have explored 
Skillman's claim using our homogeneous models and varying U (and Rj) while keeping 
Qh fixed. We recall that, as described in §2, each of our Hii region models is defined by 
fixed values of Qh for given choices of the gas density and filling factor and the ionization 
parameter (U) is defined as the ratio of the ionizing photon density to the gas density, 
evaluated at the inner radius. We emphasize that the Qh range covered by our models is the 
same as that of Skillman's models. In all cases we find that the models with fixed Qh but 
different U yield rj values which differ by less than 10% and icfs which differ by less than 
0.02%. Exceptions do occur if the U value is too low, corresponding to a very large value 
of Rj. These latter cases correspond to unrealistic models in that the geometrical width 
of the ionized region is much smaller than the inner radius (R^— Rj ^ Rj). Nevertheless, 
even for these cases (low U, fixed Qh) the icf differs from that of the higher U models by 
less than 0.5%. Thus, if U is varied while keeping Qh (and the density) fixed, the radiation 
softness parameter rj is virtually unchanged, and there is a well-defined relation between the 
ionization correction factor and the Vilchez & Pagel (1988) radiation softness parameter. 



3.2. Homogeneous Models 



First consider our results for the homogeneous models (the sohd lines in Fig. 1). For 
our fiducial models we fixed the density at njj = 10 cm""^ and varied from 7.5 x 10^^ 

to 7.5 X 10^^ s~^. The lower bound on is set by the requirement that the Hii region 
be bright enough to be observable. The arrows on the sohd lines indicate the effect of 
increasing Qh- By varying our choice of the Qh and the Hii region density we found that 
the solid (and dashed) lines arc actually the loci of constant values of the combination 
n^Q,H- In addition, although our fiducial models assume a filling factor of unity (e = 1), we 
experimented with several choices of fiUing factor (e < 1) and verified that these models 
also he along the fines in Figure 1, with a decrease in e corresponding to an increase in 
n^Q,H- As noted in §3.1 above, the range in Qh covered by our models explores the same 
ionization parameter range that Skillman (1989) covered in his investigation. 

As anticipated for the case where the nebula is ionized by a young starburst with 
hot stars providing a hard radiation spectrum, the H ii region models predict a "reverse" 
ionization correction {icf < 1); i.e., neutral hydrogen is present where the helium is still 
ionized. As n^Qn increases (and/or the filling factor decreases), r] decreases (from 77 ^ 5 or 
log 77 ^ 0.7) and the He"*" and H+ regions more nearly coincide {icf — 1) with the icf - rj 
relation being traced out by the solid fine in the lower, left-hand corner of Figure 1. 

In contrast to the young starburst case, the softer spectrum of the ionizing radiation 
present 2.5 Myr after the initial burst leads to very different behavior (see the upper, 
right-hand corner of Fig. 1). In these cases neutral helium is present in the hydrogen ionized 
zone {icf > 1). At first, as Qh increases r] barely changes, decreasing only very slightly, 
while the icf decreases noticeably (in contrast to the young starburst case). However, after 
turning the corner at the "elbow" located at log rj ^ 0.8, 77 then increases with increasing 
Qh while, now the icf remains relatively constant. 
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These distinct behaviors of H ii regions ionized by hard and by soft spectra were also 
found by Pena (1986), by Dinerstein & Shields (1986) and by PSTE. Note in particular that 
for regions ionized by the soft spectra of "old" starbursts the radiation softness parameter, 
77, is bounded from below (log 77 ^ 0.8). Since Hii regions with such large values of 77 
are normally excluded from analyses whose goal is the determination of the primordial 
abundance of helium, such old starbursts cannot dominate the ionizing radiation of the 
regions employed in such analyses. Nonetheless, as noted by Dinerstein and Shields (1986), 
some contamination from regions ionized by such soft spectra could effect the overall 
ionization correction. We explore this possibility below in §3.4. 



3.3. Inhomogeneous Models 

Now consider the non-composite inhomogeneous models (C/V) whose construction 
was described in §2. Since the results of the inhomogeneous "valley" models are virtually 
indistinguishable from those of the corresponding homogeneous models, we have not shown 
them in Figure 1. For the "condensation" models we found that the greatest deviation 
from the corresponding homogeneous case occurs when the condensation is located at a 
radius equivalent to an optical depth of r = 0.04. At this distance the radiation field is 
still sufficiently strong to have been only partially absorbed by the high density gas and the 
regions shadowed by the condensation arc still partially ionized, with the amounts of H° 
(for the high temperature cluster) and He° (for the low temperature cluster) determining 
the icf value. It is these cases (r = 0.04) which are shown by the dashed hues in Figure 1. 
Thus, all of our models, the non-composite C/V models as well as the composite models 
described below, lie between the homogeneous cases (solid lines) and these "extreme" 
inhomogeneous models (dashed hues). 

For our Monte Carlos we used the H/V/C models as building blocks in the construction 
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of a suite of composite inhomogeneous models. First, we fix the number of ionizing photons, 
Qh: and the density, n, and archive the results of a set of corresponding H/V/C models 
with the valleys and condensations located at the different optical depths (see §2). A 
single composite inhomogeneous model is created by assigning a weight proportional to the 
solid angle, Qj, occupied by each H/V/C model and summing the weighted contributions 
of each model. The Qi are chosen randomly, subject to the constraint that their sum is Att. 
This procedure is then repeated 10,000 times and the Monte Carlos are run for different 
choices of Qh and n. 

Not surprisingly, the results of our Monte Carlos (for the young starburst case) all lie 
in the shaded region in Figure 1, between the icf — rj relation for the homogeneous models 
(solid curve) and that for the "extreme" condensation (r = 0.04) models (dashed curve). 
Although the behavior for the case of the old starburst is more complicated, it is clear that 
these composite inhomogeneous models will always have "large" values for the radiation 
softness parameter and cannot, by themselves, describe the Hii regions selected for probing 
the primordial helium abundance. 



3.4. Contamination By "Old" Hii Regions 

As we have seen, H ii regions ionized by young starbursts with hard spectra always have 
a "reverse" ionization correction {icf < 1) and are limited to relatively low values of the 
radiation softness parameter (log rj ^ 0.7). In contrast, regions ionized by a softer spectrum 
will have icf > 1, but also log r] ^ 0.8. Although the latter would normally be eliminated 
from consideration by an r] cut, it is possible that a contribution from such regions could 
contaminate the helium abundances derived from observations (Dinerstein & Shields 1986). 
To test this possibility it would be necessary to construct models of superposed Hii regions 
for each observed H ii region and to constrain the parameters of the mixture by demanding 
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that these composite models reproduce all the observed line strengths. As an illustrative 
example, Dinerstein and Shields (1986) have done something simpler. For NGC 4861 they 
superposed a "hot" region (T* = 55,000 K) and a "cool" region (T* = 35,000 K), fixing 
the ionization parameter for each region by the arbitrary requirement that each component 
reproduce the observed I(A6300)/I(H/3) ratio. The relative weights of each component were 
then fixed by the requirement that the resulting I(A3727)/I(H,5) ratio have the observed 
value. In this case they found that the low-T* component contributes 12% of the observed 
H/3 luminosity. For this composite model they found icf — 1.09 (and log 77 = 0.39) in 
contrast to the value icf = 0.99 (and log t] = 0.47) for their simple, non-composite model, 
and they noted that ''low observed values o/He"'"/H"'" cannot be interpreted as low He/H 
with complete confidence" . However, it should be remarked that neither of their models is 
consistent with all their observed line ratios. For example, neither their one-component 
nor their two-component model reproduces the observed value of the radiation softness 
parameter: log rjobs — —0.07. Notice from Figure 1 that for our one-component homogeneous 
and inhomogeneous models, the icf ^ 0.99 for log rj —0.07. 

As suggested by the failure of Dinerstein and Shields (1986) to account for all the 
observed line ratios, it may not be so easy to confuse a composite "young/old" Hii region 
with a single (albeit inhomogeneous) Hii region. Various line ratios can be the key to 
distinguishing between the two possibilities. Using our young (t = 0) and old (t = 2.5 Myr) 
H II region models we have constructed a suite of composite models whose results are easy to 
understand with reference to Figure 1. For almost all combinations one or the other region 
will dominate and the icf — rj results will lie close to those of the dominant case (i.e., the 
one with the largest Qij). Given the separation in log 77 between the young and old regions, 
it is unlikely that an observed H 11 region will be confused by a composite region dominated 
by a component with the "wrong" value of log 77. Since virtually all of the Hii regions 
selected for determination of the primordial helium abundance have low values of log 77 
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(see the data utilized in the next section), this suggests that most of the viable composite 
models will have icfs similar to those for our t = case {i.e., icf ^ 1). Only when the 
number of ionizing photons (Q//) for each of the two regions (t = and t = 2.5 Myr) are 
comparable may it be possible for a composite region to masquerade as a single region. We 
have therefore concentrated on such "democratic" composite models and found that they 
occupy a region in the icf - log r] plane with log r] ^ 0.4 and icf ^ 1.04. For the IT data 
which we use in our quantitative analysis in the next section, log r] ^ 0.37 and interloping 
regions are unlikely to have played a major role. Nonetheless, we have pursued the question 
of whether such composite regions could contaminate the icf inferred for H ii regions which 
are observed to have log r] ~ 0.4. 

As a first step we have compared the results of the t = 0, Qh = 7.5 x 10^''' model 
(n = 10 cm~^ and e = 1) for which log r] = 0.37 and icf = 0.956 with those for a composite 
model of two regions (t = and t = 2.5 Myr) each with Qh — 7.5 x 10^^ {n — 10 cm~^, 
e = 1) for which log rj — 0.42 and icf — 1.036. Notice that if e < 1, similar r) and icf values 
are obtained with higher n^Qn- Although the rjs of the single and composite regions are 
very similar, we find that the individual ion ratios 0"'""'"/0"'" and S"'""'"/S"'" differ by factors of 
6-7 between the simple and the composite models. To explore this further, we compared 
the ratios of the intensities of various emission lines to that of H/3 for the simple model (t 
— 0, Qh = 7.5 X 10'^'' s~^) with those for a series of composite models with Qh varying 
from 7.5 x 10^^ to 7.5 x 10^^ s~^. All of these composite models have log rj ^ 0.4 and 
icf 1.04. In the composite models the [Oii](A3727) to H/3 ratio varies from 18% to 66% 
of that in the simple model, the [N ii](A6584+6548) to H/3 ratio is 16% to 60% of that in 
the simple model, and the [Sii](A6717+6731) to H/3 ratio is 3% to 41% of that in the simple 
model. In contrast, the [Oiii](A5007+4959) to H/5 line ratios in the composite models 
exceed those in the simple model by factors of 1.4 to 2.3. Prom our analysis it seems clear 
that these crucial line ratios can be key to unmasking composite H ii regions masquerading 
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as simple Hii regions (see, also, Pena 1986). 

4. Corrections To The Helium Abundance For Observed Hii Regions 

As we have done in our previous analysis of the effect of temperature inhomogeneities 
on the primordial helium abundance inferred from the observational data (Steigman, Viegas 
& Gruenwald 1997), we quantify the effect on Yp of an icf which differs from unity using a 
sample of H ii regions from the literature. The data of Izotov and Thuan (IT) is especially 
well-suited to our task. For 41 of the 45 Hii regions presented by IT, they provide helium 
abundances and data from which we may estimate the radiation softness parameter r]. 
Although one of their Hii regions, UM 311, is not really in the category of "metal-poor", 
we have verified that whether or not we include it in our analysis has no significant effect 
on the difference in the derived value of Yp using our calculated icf as compared to the 
IT choice icf — 1 for all regions. For the IT sample, —0.25 < log r] < 0.37 suggesting 
(see Fig. 1) an "average" icf ^ 0.96 ± 0.02 and a potentially large reduction in Yp, 
AYp f» —0.0072 ± 0.0036. To quantify this reduction we have run a series of Monte Carlos 
using the IT data set and our computed icf — r] relations shown in Figure 1. 

For each Hii region in the IT sample we have 77, the oxygen abundance 0/H, and Y 
for icf = 1. For this fiducial data set we find the best linear fit to Y versus 0/H in order 
to derive Yp. We have also experimented with using averages of Y for the 15-20 lowest 
oxygen abundance regions to bound Yp and find the difference in our derived AYp to be 
negligible compared to that found using the linear fit. We emphasize that here we are 
not interested in the actual value of Yp per se. Rather, we want to find the reduction 
in Yp due to icf < 1. In the Monte Carlos, for each of the observed regions we use the 
value of log rj from the data to randomly choose an icf between the minimum value for 
a condensation located at r = 0.04 ("extreme" inhomogeneous case) and the icf for the 
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homogeneous case. In choosing between these two extremes we have experimented with 
three probabihty distributions: equal probabihty for the icf to he between the two extremes 
and decreasing/increasing probabihty for the extreme inhomogeneity. For each of the (41) 
IT regions we correct their value of He+/H+ with the randomly chosen icf and compute 
the ic/-corrected He+/H+ ratio which we use to find the ic/-corrected Y. This new set 
of (41) Y, 0/H pairs is used to find the corresponding Yp from the linear Y versus 0/H 
fit. For each choice of icf probability distribution we redo this 10,000 times. Our results 
appear as histograms for the distributions of AYp in Figure 2. It is clear from Figure 2 
that our principal conclusion, —AYp fa 0.003 — 0.004, does not depend on the choice of icf 
probability distribution. Our Monte Carlos thus suggest that by ignoring the correction for 
He - H ionization, IT have overestimated the primordial helium mass fraction by 0.0034 
± 0.0003. For their sample, IT derive (see their Table 7) Yp (IT) = 0.2443 ± 0.0015 (or, 
0.244 ± 0.002). Including our estimate of the true icf we suggest their value of Yp should 
be reduced to Yp(VGS) — 0.241 ± 0.002. We note that icf = 1 has also been assumed for 
the Hii regions in the data sets used by OS and OSS, so their estimates of Yp should likely 
be reduced by a similar amount to that which we have found for the IT sample. 

5. Conclusions 

Relatively accurate helium abundance determinations are currently available for nearly 
100 low-metallicity extragalactic Hii regions (see, e.g., OS; OSS; IT). These data enable 
estimates of the primordial abundance of helium to unprecedented statistical accuracy. For 
example, IT quote AYp = ± 0.0015 (or, ± 0.002) for the regions they have observed and 
analyzed. At this level of uncertainty unanticipated systematic errors may overwhelm the 
statistical uncertainties (see, e.g., Davidson & Kinman 1985; PSTE; Skillman et al. 1994; 
ITL; Peimbert 1996; Steigman, Viegas & Gruenwald 1997; Skillman, Terlevich & Terlevich 
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1998; IT). In this paper we have revisited the question of the hehum ionization correction 
and its contribution to estimates of Yp. Virtually all the low-metallicity Hii regions 
employed in the quest for Yp have "hard" spectra (low r]) and, almost always the ionization 
correction has been ignored or, rather, assumed to be unity (no correction). We have 
seen, as have others before us {e.g., Pena 1986, Dinerstein & Shields 1986; PSTE), that 
there is a reverse ionization correction {icf < 1) for regions ionized by such hard spectra. 
Inhomogeneities, surely present in these regions, only serve to exacerbate this correction 
{i.e., to reduce the icf). In our comparison with real data (IT), we suggest that their 
estimate of Yp = 0.244 ± 0.002, should be reduced to Yp = 0.241 ± 0.002, a correction 
larger than the quoted statistical error. Of course, this latter estimate still ignores other 
possible sources of systematic error such as temperature fluctuations (Peimbert 1971; 
Steigman, Viegas & Gruenwald 1997) and underlying stellar absorption (see ITL; Skillman, 
Terlevich & Terlevich 1998; IT). 

In the current precision era of cosmology predictions and observations are achieving 
unprecedented levels of statistical accuracy. This is certainly the case for primordial 
nucleosynthesis in the standard (isotropic, homogeneous, three flavors of light neutrinos, 
etc.) hot big bang cosmological model (BBN). For example, at a fixed value of the baryon 
density (or, equivalently, the baryon-to-photon ratio) the uncertainty in the predicted 
primordial abundance of helium is at the level of ± 0.0005 (Hata et al. 1996; Buries 
et al. 1999; Olive, Steigman & Walker 1999; Lopez & Turner 1999). Furthermore, the 
variation of Yp(BBN) with baryon density is very slow so that if deuterium is used as 
a baryometer to constrain the baryon density, AYp(BBN) f» 0.006(Ay2/y2) where y2 = 
(D/H)p. The uncertainty in the BBN-predicted primordial deuterium abundance is of order 
8% (Hata et al. 1996) or less (Buries et al. 1999), a level comparable to that claimed for 
the uncertainty in the value inferred from observations of deuterium along the lines-of-sight 
to two high-redshift, low-metallicity (hence, very nearly primordial) QSO absorption line 
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systems (Buries & lytler 1998a,b): y2{obs) = 3.4 ±0.25 x 10"^ The BBN helium abundance 
which corresponds to this deuterium abundance is Yp(BBN) = 0.247 ± 0.001. Notice that 
while this expected abundance is consistent with the IT determined value (0.244 ± 0.002) 
it is higher than our IT-corrected estimate (0.241 ± 0.002) by some 3a. 

Strictly speaking, our quantitative estimate of the reduction in Yp associated with the 
ionization correction factor apphes only to the data of IT. Nonetheless, it is anticipated 
that Yp inferred from the data assembled by OS should be similarly reduced. Since the OS 
estimate of Yp (excluding the NW region of IZwlS which is suspected of being contaminated 
by underlying stellar absorption) is even lower than that of IT (Yp(OS) = 0.234 ± 0.003) 
the discrepancy between theory and data is even larger. Of course it could be that the 
Buries and Tytler (1998a,b) estimate of the primordial deuterium abundance is too low 
(see, e.g., Webb et al. 1997; Levshakov, Kegel and Takahara 1998a,b, 1999). Nonetheless, 
this example provides an object lesson highlighting the importance of careful estimates of 
systematic errors in the current era of precision cosmology. 

In deriving the helium abundance from observations of Hii regions it is necessary 
to correct the emission-line data for various effects {e.g., coUisional excitation). In this 
paper we have explored the correction for unseen hydrogen and/or helium using detailed 
photoionization models of H ii regions ionized by realistic radiation spectra from different 
(young/hot and old/cool) stellar clusters. For regions ionized by young, hot stars the 
radiation softness parameter is small (log t] ^0.7) and there is a systematic, reverse 
ionization correction {icf < 1). In contrast, for regions ionized by an older stellar population 
with fewer hot stars the radiation softness parameter is large (log 77 ^ 0.8) and icf > 1. 
Thus, the correction for ionization is systematically correlated with the magnitude of the 
radiation softness parameter. Of course, if this correction were very small it would have a 
negligible effect on the helium abundance determination and could be neglected. Our results 
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suggest this is not the case. For the H ii regions in the IT data set there are reductions 
in individual Y values which could be as large as 0.01 or even larger. These reductions, 
which are comparable to (or, even larger than) typical statistical errors in Y for individual 
Hii regions (see, e.g., OS, OSS, IT), should be compared to the theoretical (SBBN) 
uncertainty which is more than one order of magnitude smaller(e.^., Olive, Steigman 
and Walker 1999). For the IT data set we find the reduction in the inferred primordial 
mass fraction is comparable to (or, even somewhat larger than) the quoted statistical 
error in Yp. The importance of including the ionization correction may be illustrated by 
the corresponding bound on N,^, the "effective" number of light neutrino species, which 
provides a measure of any deviation from the standard model energy density at the epoch 
of BBN (Steigman, Schramm and Gunn 1977). For SBBN, = 3, while the presence 
of "new" , light particles (beyond the standard model of particle physics) would permit 

> 3. In a recent analysis utilizing the IT determination of Yp and the Buries and Tytler 
(1998a,b) deuterium. Buries et al. (1999) find a 2cr upper bound of N,^ < 3.2 (see, also, 
Olive, Steigman and Walker 1999). Incorporating the reduction in Yp from the ionization 
correction will reduce this bound to Nj, < 2.9, posing a potential challenge to SBBN. While 
there undoubtedly remain other, as yet unquantified, systematic errors whose magnitude 
may be larger than the ionization correction, it seems clear that the systematic reduction 
for ionization cannot be ignored. 
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Figure Captions 

• Figure 1: The icf— log rj relations for Hii regions ionized by "young" (t = 0) and 
"old" (t = 2.5 Myr) star clusters. The results for the young star clusters occupy the 
lower, left-hand part of the figure, while those for the older star clusters are found 
in the upper, right-hand part of the figure. The solid lines are for the homogeneous 
models, while the dashed lines are for the "extreme" inhomogeneous models (see the 
text). The shaded region between the sohd and dashed fines for the t = case are the 
results of Monte Carlos for composite inhomogeneous models. 

• Figure 2: Distribution of the offsets in the primordial hefium mass fraction, Yp, 
inferred from Monte Carlos using tfie IT data and our icf — ?] results for tfiree 
different probability distributions whicfi interpolate between the homogeneous and 
"extreme" inhomogeneous models (see the text) . Prom top to bottom the probabifity 
distributions favor the homogeneous models, are neutral, and favor the inhomogeneous 
models. 
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